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ABSTRACT 


The  Wavenumber  Technique  (WI)  is  a  relatively  new  method 
of  underwater  sound  transmission  analysis.  One  aspect, 
source  depth  determination,  is  studied  to  evaluate  its 
validity  and  test  environmental  and  acoustic  sensitivity. 
The  horizontal  wavenumber  spectrum  is  analyzed  to  determine 
null  spacings  in  wavenumber  space,  which  indicates  source 
depth  by  the  Lloyd's  Mirror  interference  affect.  Comparison 
of  this  theory  with  cases  of  an  isospaed  sound  profile, 
fully  absorbing  bottom,  and  flat  totally-reflecting  surface 
shows  excellent  agreement  for  several  parametric  variations. 
Cases  with  realistic  sound  speed  profiles  and  partially 
absorbing  bottoms  generally  agree  with  theory,  but  a  dis- 
tinct bias  is  observed.  Source  deprh  determination  curves, 
which  relate  the  scaled  wavenumber  specrral  intensity  null 
spacing  to  the  sourc?  depth,  are  presented  for  comparison 
with  theory  and  sensitivity  analysis.  An  example  is  given 
for  suggested  application  of  source  depth  determination. 


TABLE   OF    CONTENTS 

I.  INTRODUCTION 17 

II.  'WAVENOcMBES    TECHNIQUE 19 

&.       GENERAL    DESCRIPTION 19 

B.  PHYSICAL    DESCRIPTION 23 

C.  WAVENUMBER    TECHNIQUE    APPLICATION       28 

D.  TECHNIQUE    LIMITATIONS     30 

III.  PROPAGATION    MODELS 31 

A.        INTRODUCTION 31 

3.        PARABOLIC    EQUATION 33 

1.  Introduction 33 

2.  Parabolic   Approxiiation 34 

3.  Brock   Algorithm 35 

a.  Description 35 

b.  Implementation 37 

U.      Subsequent    Modifications U1 

a.  Introduction U1 

b.  Tatro    Modifications U1 

c.  Stieglitz   at    al.  Modifications U2 

5.  SSFFT   Evaluations    45 

6.  Environmental  Sensitivi-y   of    the   PS   ....  49 
C.        FAST    FIELD     PROGRAM 52 

1.      General   Description 52 

5 


2.      Comparison    of  FFP    and  PE SU 

D.  NORMAL    MODE    THEORY 55 

1.  General   Description   55 

2.  Comparison    of  NM   and  PE 57 

E.  FINITE    DIFFERENCE   METHODS 58 

1.  General   Descrip-ions 58 

2.  Comparison    of  FD   and  DDE    with    PE 60 

IV.  ANALYSIS    OF    MT 62 

A.  GENERAL    DESCRIPTION    62 

B.  ISOSPEED    CASES 87 

1.  Water  Column   Dsp-h    Variation      87 

2.  Range   Variation 89 

3.  Frequency   Variation 89 

U.      Receiver  Depth   Variation 96 

5.       Analysis  of   Isospeed  Cases      96 

C.  REALISTIC    PROFILE   CASES    100 

1.  Spatial   Variation    100 

2.  Temporal  Variation      101 

3.  Range    Variation 105 

U.      Receiver  Depth   Variation 105 

D.  COMPOSITE    COMPARISON 105 

?.  CONCLOSIONS 111 


LIST  OF  REFERENCES 115 

INITIAL  DISTRIBUTION  LIST   118 


LIST   OF    TABLES 

TABLE   I.  Scenario    Specifications      64 

TABLE  II.  Location   Informa-ion   for    Realistic   Profiles  68 

TABLE  III.        Coipariscn   Specifications 86 


8 


LIST   OF    FIGURES 

Figure    1.         Wavenumber  Techniqas    Flow    Diagram 21 

Figure  2.         Lloyd's   Mirror  Effsct   Geometry 26 

Figure  3.         WT   Application   for  Source   Depth 

De-ermination 29 

Figure  U,         Space   and  time  relationships   of   actual 

profiles 67 

Figure   5.         Isosoeed    Profile   Waven umber   Spectrum    for    SD 

of    200    ft 70 

Figure  6.         Isospeed    Profile   Wavenumber   Saectrum    for    SD 

of    300    ft " 71 

Figure  7.        Isosoeed    Profile   Wavenumber  Spectrum   for   SD 

of    500    ft 72 

Figure  8.         Isosoeed    Profile   Wavenumber  Spectrum   for   SD 

of    800   ft 73 

Fiaure  9.         Isosoeed    Profile   Wavenumber   Soectrum    for    SD 

of    1500    ft ■ 7U 

Fiaure   10.      Isosoeed    Profile    Wavenumber   Spectrum    for    SD 

of  3000  ft 75 

Figure  11.   Realistic  Profile  Wavenumber  Spectrum  for  SD 

of    200    ft 76 

Figure   12.      Realistic   Profile   Wavenumber   Spectrum    for    SD 

of  300  ft 77 

Figure  13.   Realistic  Profile  Wavenumber  Spectrum  for  SD 

of  500  ft 78 

Figure  l^i.   Realistic  Profile  Wavenumber  Soectrum  for  SD 

of  800  ft : 79 

Fiaure  15.   Realistic  Profile  Wavenumber  Spectrum  for  SD 

of    1500    ft 80 

Figure    16.      Realistic  Profile  Wavenumber  Sp=ctrum    for    SD 

of   3000    ft 81 


Figure  17.   Isospeed  Profile  Wavenumber  Spectrum  for  SD 

cf  50  ft ^ 82 

Figure  18.   Example  of  Source  Depth  Detsrmination  Curve.   85 

Figure  19.   Isospeed  Profile  SD  Detarmination  Curve:  Z 

Variation  -  1 88 

Figure  20.   Isospeed  Profile  SD  Determination  Curve:  Z 

Variation  -  2 90 

Figure  21.   Isospeed  Profile  SD  Determination  Curve:  Z 

Variation  -  3 91 

Figure  22.   Isospeed  Profile  SD  De-armina*:  ion  Curve:  Z 

Variation  -  U 92 

Figure  23.   Isospeed  Profile  SD  Detsrminat ion  Curve:  R 

Variation  -  1 93 

Figure  24.   Isospeed  Profile  SD  Determination  Curve:  R 

Variation  -  2 9a 

Figure  25.   Isospeed  Profile  SD  Detsrmination  Curve:  R 

Variation  -  3 95 

Figure  26.   Isospeed  Profile  SD  Det=rminat ion  Curve:  F 

Variation  -1 97 

Figure  27.   Isospeed  Profile  SD  Determination  Curve:  F 

Variation  -2 98 

Fiaure  28.   Isospeed  Profile  SD  Determination  Curv=:  RD 

Variation 99 

Figure  29.   Spatial  Variation  (^,B):  SD  Determination  and 

AXBT  Curves 102 

Figure  30.   Spatial  Variation  (C,D):  SD  0= termination  and 

AXBT  Curves 103 

Figure  3  1.   Temporal  Variation  (A,C):  SD  Determination 

and  AXBT  Curves 10U 

Fiaure  32.   Temporal  Variation  (B,D) :  SD  Determination 

and  AXBT  Curves 107 

Figure  33.   Range  (R)  Variation  for  Profile  A 108 

Figure  3U.   Receiver  Depth  (RD)  Variation  for  Profile  A.  109 


10 


Figure  35. 


Conposite  Comparison 
Realistic  Cases.  . 


for  Isospeed  and 


110 


11 


TABLE   OF    SYMBOLS 

a  =  attenuation   coefficient 

3  *  scaled   wavenumber 

1^;  =  envelope   function 

\  =  reference    wavelength 

0 

ui  =  angular    frequency 

p  =  density 

9  =  beam   elevation     (depression)    angle 

5  =  central-diff ere  rx:e   operator 

^  =  3.  1U1593 

7  =  Laplacian   operator 

3  =  partial   derivative   operator 

Z  =  summation    operator 

/  =  integral  operator 

v''~  =  square   root    operator 

<<  =  much    less   than 

C  =  sound   speed 

C  =  minimum    sound  sceed 

0 

C  =  normal    mode   phase    soeed 
P 

Z  =  Lloyd's    Mirror    pressure   field  in   raage   space 

f  =  frequency 

?  =  Lloyd's    Mirror    pressure   field   in   wavenumber   space 

G  =  Green's    function 

12 


H  '  Hankel   function 

0 

I  =  spectral   intensity 

Im^,  *  imaginary    pressure    transform 

J^  =  zero   order  Bess  el    function 

K  =  total   wavenumber 

k  =  horizontal   wavenumber  component 
r 

k  =  initial    spectrum   wavenumber 

k  =  vertical    wavenumber   component 

Ak  =  wavenumber   increment 

K   ^,=  maximum    wavenumber 

max 

M  =  complex    modified  index   of    refraction 

„  =  index  of   refraction 

n  =  modified   index    of   refraction 

n,  =  water/bottom   interface   index  of  refraction 

D 

N  =  spectrum   index 

P  =  acoustic   pressure 

:?  =  ra  n  g  e 

AR  =  range  increment 

^e^  =  real  pressure  transform 

u  =  wave  speed 

Z  =  wa-er   column   depth 

Z  =  modified    water    column  depth 

AZ  =  dep-h   increment 


13 


Z,  =  receiver  depth 

Ze  =  source  depth 

2  -   maximum  depth  in  transform 
max 


1U 


ACKNOWLEDGEMENT 

Dr.  A.  E.  Coppens  showed  great  kindness  by  accepting 
supervision  of  this  research  when  time  was  short.  His  con- 
cern, understanding  and  direction  led  to  an  investigation 
which  is  believed  to  te  successful.  He  can  both  profess  and 
teach.  This  research  and  my  personal  development  benefited 
greatly    from   this    unique    ability. 

LCDR  C.  E.  Dunlap,  USN  (Ret)  was  very  supportive  and 
encouraging  during  my  entire  thesis  research  effort.  I 
appreciate  his  assistance  with  this  project  and  his  excel- 
lent work  with  the  NPS  Environmental  Acoustics  Research 
Group     (EARG)  . 

Richard  Lauer  of  the  Naval  Ocean  Rssearch  and  Develop- 
ment Activity  (NORDA)  developed  the  MT  coicept  which  was  the 
basis  for  this  research.  My  -hanks  to  him  for  allowing  -his 
investigation  of  his  idea  and  his  assis'rance  during  the 
?-udy.  Tom  Lawrence  of  Ocean  Data  Systems,  Inc.  (ODSI)  at 
NORDA  was  my  principal  contact  during  this  research  and  fur- 
nished i-s  foundation.  I  appreciate  the  iiany  hours  of  tele- 
phone conversation  and  his  liaison  at  NORDA.  Richard  Evans, 
also  of  ODSI  and  NORDA,  provided  great  assistance  in  devel- 
oping the  computer  software  and  my  understanding  of  the  mod- 
els. 

15 


I  also  thank  C.W.  Spofford  and  Eleaa^r  Holmes  of  Science 
Applications,  Inc.,  John  Locklin  and  Gil  Jacobs  of  ODSI,  and 
Jim   Davis   of   Planning   Systems,    Inc.    for   their   assistance. 

The  personnel  of  the  NPS  Dudley  KnDx  Library  Research 
Reporxs  Division  and  the  W.R.  Church  Computer  Center  were 
also   very   helpful   during    this   research. 

My  -hanks  also  to  my  family  for  thsir  patieice  and 
understanding,  and  to  my  colleagues  in  the  EARG  for  their 
advice  and   assistance. 


16 


I.        INTRODUCTION 

The  Wavenumber  Technique  (WT)  is  a  new  method  for  the 
solution  of  several  problems  related  to  the  analysis  of 
underwater  sound  transmission  (Lauer ,  1979)  .  Its  basis  is  the 
analysis  of  acoustic  propagation  in  the  wavenumber  domain, 
which  has  been  described  by  DiNapoli  (1971)  and  was  an 
in-ermediate  step  in  the  computation  of  transmission  loss 
versus  range  in  the  Fast  Field  Progran  (FFP)  .  DiNapoli 
(1977)  later  used  the  wavenumber  domain  to  study  the  impe- 
dance of  the  ocean-bottom  interface.  Lauer  (1979)  first 
introduced  the  WT  and  its  applications  to  passive  localiza- 
tion and  tracking  and  multipath  decomposition.  He  developed 
a  range  prediction  curve  based  upon  peak  wavenumbers  in  the 
transform  spectrum  for  successive  range  increments  from  an 
entire  pressure    range   field. 

These  investigations  presented  highly  beneficial  and 
promising  results  from  the  wavenumber  depiction,  which  indi- 
cates that  additional  research  is  appropriate.  One  particu- 
lar aspect  of  the  WT  -  determination  of  source  depth  -  will 
be  investigated  here  in  order  to  evaluate  its  validity  in 
comparison    with      the    Lloyd's    Mirror   effect.  Its   relation- 

ships to  oceanic  and  acoustic  variations  will  also  be 
examined . 

17 


The  wavenumber  technique  as  first  described  by  Lauer 
(1979)  was  based  upon  input  from  the  FFP,  but  it  has  been 
adapted  for  use  with  the  Parabolic  Equation  (PE)  method  of 
propagation  loss  determination.  The  PE-based  WT  is  currently 
being  investigated  by  Lauer  and  others  at  the  Naval  Ocean 
Research  and  Development  Activity  (NORDA)  and  was  the  model 
used  in  this  research  at  the  Naval  Postgraduate  School 
(NPS) .  Descriptions  of  the  PE  and  WT  will  be  presented  to 
display  their  sensitivities  to  variations  of  several  model 
parameters.  Also,  expected  relationships  of  the  WT  to  other 
current  propagation  prediction  models  will  be  discussed  to 
consider  how  results  might  vary  based  upon  different  inputs. 
Particular  attention  will  be  given  to  algorithm  differences, 
degrees  of  approximation,  and  treatment  of  the  ocean  bottom. 

The  PE-based  WT  will  then  be  exercised  for  several 
parameter  variations  to  evaluate  qualitatively  the  consis- 
tency cf  the  technique.  This  study  will  also  demonstrate 
the  response  of  the  PE-based  WT  to  geometric  and  actual 
oceanic  conditions.  Prediction  curves  for  source  depth 
d= terainat ion  based  upon  these  variations  will  be  presented 
and  discussed. 
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II.       WAVENUMBSR    TECHNIQUE 

A.        GENERAL    DESCRIPTION 

The  WT  is  a  natural  alternative  for  signal  depiction 
because  it  provides  information  on  the  directions  of  energy 
arrivals.  It  allows  reasonably  direct  physical  interpreta- 
tion of  ocean  acoustic  processes  and  their  relationships  to 
environaental  conditions  (Lauer, 1979) .  Lauer  further  states 
that  the  wavenumber  spectrum  plot  is  clearer  and  more 
informative  than  typical  curves  of  propagation  loss  versus 
range.  The  wavenumber  K  relates  the  angular  frequency  and 
the    sound   speed   by 


K     =      00      =      2TTf 

—       J 

c        : 


and  it  has  horizontal  and  vertical  components 


K  =  k   +  k   . 
r  z 


The    reference      wavenuirber    is   determined      by    the      sound   speed 


minimum 


K       =    CO 
o         — 

C 

0 
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The  WT  will  first  be  described  withoat  reference  to  any 
particular  propagation  model  to  present  the  general  nature 
of  the  algorithm.  Fig.  1  is  the  flow  diagram  for  the  WT. 
The  left  and  right  sides  of  this  flow  diagram  correspond  to 
the  intensity  and  waven umber  axes,  respectively,  of  the 
spectrum  plot  that  will  result  from  the  WT.  This  figure  and 
the  following  discussion  were  developed  from  Lauer  (1979  and 
1982a)  and  augmented  by  several  telephone  conversations  with 
R.  Lauer,  T.  Lawrence,  and  R.  Evans  of  NORDA  between  July 
and    December,    1982. 

The  WT  calculates  spectral  intensity  from  the  real  and 
imaginary    parts   of  .the  acoustic    pressure 


?    =    ^(R,Z)    e^^^ 


The  complex  psi  versus  range  field  at  a  specified  depth  is 
the  propagat.icn  model  product  which  detariines  the  pressure. 
The  complex  pressure  field  can  be  modifiad  for  effects  such 
as  attenua-cicn,  depending  upon  whether  ii=alized  or  realis- 
tic conditions  are  being  investigated.  A  Fast  Fourier 
Transform  (FFT)  is  then  applied  to  yield  transformed  complex 
pressure.           The   spectral      intensity    is      ds-ermined   at      each 
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wavenumber   increment      from  the  transformed      complex   pressure 
field   by 


I^    -   Re^2    +    Iin^2    . 


This  intensity  field  is  scaled  arbitrarily  for  satisfactory 
display  of  the  relative  maxiia  and  minima  to  define  the 
intensity  nulls. 

Additional  output  parameters  from  the  PE  which  are 
required  for  the  WT  include  tha  range  stap,  the  number  of 
range  points,  frequency,  minimum  sound  speed  and  field 
dsp-h.  These  are  the  selected  variables  noted  in  figure  1. 
The  range  step  is  selected  to  provide  sufficient  resolution 
*o  determine  the  wavenumber  null  spacing  on  the  spectrum 
pio-s.  The  initial  radial  coaponent  of  the  wavenumber  for 
the  spectrum  is  calculated  from  DiNapoli  (197j1)  by 

AR 


where 


K   =  k    +  k 
0     r  0     z  0 
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As  will  be  seen  in  section  II-B,  if  a  "scaled"  wavenum- 
ber  beta  is  used  for  the  horizontal  axis,  analysis  is  sim- 
plified since  the  nulls  in  the  spectra  should  be  more  evenly 
spaced  (Lauer, 1979)  .  The  physical  meaning  of  beta  will  be 
described  in  the  next  section.  The  resulting  scaled  wave- 
number  spectra  are  then  plotted.  The  wavenumber  null  spacing 
for  each  incremental  source  depth  is  determined  from  each 
scaled  spectrum.  The  source  depth  determination  curve  is 
generated  by  plotting  these  null  spacings  as  a  function  of 
source  depth,  as  suggested  by  Lauer  (1979).  The  null  spac- 
ing cf  an  analyzed  received  signal  could  be  compared  with 
this  determination,  curve  to  infer  the  source  depth  of  the 
signal.  This  application  will  be  discussed  later.  The 
source  depth  determination  curve  will  b=  specific  for  the 
acoustic  and  oceanic  descriptions  of  the  medium. 

E.   PHYSICAL  DESCRIPTION 

The  preceding  description  was  intended  to  provide  the 
overall  concept  of  the  «JT  wher2as  what,  follows  will  explain 
the  physical  basis  for  the  technique.  The  underlying  princi- 
ple cf  the  wT  as  described  by  Lauer  (1932a)  is  the  Lloyd's 
Mirror  affect,  which  gives  the  acoustic  field  by 

7-/-  -)  N   -  i(^t 
D  =  L  (  R )  e 
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where 


E(R)    =    exp    {iK    /  R^    +    (Z^    -    Zg)^}      -    exp    { iK    /  R^    +    (Z^^   +    Z^)^} 


/  r2    +    (Z^   -    Zg)^  /  R^    +    (Zj^    +    Zg)^ 


for  a  point  source  in  a  semi- infinite  aedium  of  constant 
sound   speed    with      a    pressure   rslease    surface.  This    is   the 

Lloyd's  Mirror  field  in  range  space,  whereas  the  Lloyd's 
Mirror  field  in  wavenumber  space,  applicable  to  this 
research,    is   given   by   Lauer      (19  82)    as 


F(K)    =    sinCSZg)    expCigZ^^)    , 


where   beta    is    defined   by 


3    =    (K    2    -   k    2)'"^ 

0  r  ^ 


and  E(R)  and  F(K)  are  related  by  the  3essel  transform  pair 


E(R)  =  /  2  F(K)  J  (KR)  K  dK 

0  0 


and 


F(K)  =  1  /  E(R)  J  (KR)  R  dR 


2  ' 


2(4 


For  this  research  the  Bessel  function  was  approximated  as  a 
trigonometric  function  so  that  the  PFT  could  be  utilized. 
This  is  considered  acceptable  since  the  source  and  receiver 
are  greatly  separated,  according  to  Lauer  (1982b),  Inspec- 
tion of  F(K)  reveals  that  the  nulls  of  the  spectrum  will  be 
equally  spaced  at  intervals  of 


This  equal  spacing  in  beta  is  the  )cey  element  of  the  appli- 
cation of  the  WT  to  source  depth  determination  (Lauer , 1 979)  . 
Only  the  direct  and  surf ace -reflected  waves  will  inter- 
act, at  a  receiver  location  in  a  simple  illustration  of  the 
Lloyd's  Mirror  effect  with  a  smooth  surface.  This  interac- 
tion will  be  destructive  or  constructive,  as  a  result  of 
their  relative  phases  (Kinsler  et  al  ,  1982)  and  will  produce 
nulls  or  peaks,  respectively,  in  the  wavenumber  spectrum. 
Fig.  2  depicts  the  Lloyd's  Mirror  effect  geometry  which  will 
be  used  to  study  tesx  case  outputs  for  null  spacing  and  sub- 
sequent source  depth  determination.  The  surface-reflected 
{S3.)  wave  travel  distance  will  increase  at  a  faster  rate 
than  the   direct  (D)   wave  distance   as  the  source   is  moved 
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Figure  2.      Llcyd's   Mirror   Effsct    Geome-ry. 
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deeper  in    the    water  cclumn    (i)       while  the   receiver   depth    (h) 
is    stationary.      This    leads  to   different   arrival    patterns   and 
nulls     at.   more      wavf=  numbers.      Thus      the      general   pattern      of 
decreased   spacing     in   beta  with      increasing    source      iepth    is 
observed . 

Eottom  interacting  waves  will  be  neglected  for  idealized 
cases.  This  will  be  accomplished  by  the  use  of  a  fully 
absorbing  bottom.  These  waves  will  be  important,  however, 
in  shallow  water  or  at  greater  ranges.  A  few  test  cases  will 
have  realistic  bottoms.  VariDus  propagation  models  handle 
the  bottom  in  different  ways,  so  the  inclusion  of  a  bottom 
will  be  discussed  later.  Cas=rS  in  shallDw  water  (less  than 
5000  feet)  or  at.  extended  ranges  (greater  than  100  miles) 
will    not    be   considered  here. 

Lauer  (1979)  observed  that  the  waveaunber  spec-ram  shows 
a  series  of  modulating  envelopes  which  enclose  spikes  that 
are  the  eigenvalues  associated  wi-^.h  the  normal  modes.  Com- 
parison of  the  wavenumber  spectrum  with  typical  transmission 
loss  curves  reveals  that  the  wavenumber  depiction  is  more 
coherently  structured.  This  leads  to  a  be'rter  physical 
understanding  of  signal  transmission  because  of  the  rela- 
tionship   between      horizontal    wavenumber    aid    ray      theory.       He 

gives  this    as 

^^    -   £   cos  9    =    2T7f    cose    . 

c  c 
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The  horizontal  wavenumber  decreases  as  theta  increases.  Thus 
the  wavenumber  axis  of  the  spectrum  plot  will  include  wave- 
numbers  associated  with  smaller  elevation  (or  depression) 
angles  in  increasing  crder  of  wavenumber  (Lauer , 1979) .  It 
is  then  expected,  for  the  case  of  a  neglected  bottom,  that 
spectrum  amplitude  peaks  could  be  observed  that  would  corre- 
spond  to   the   direct   and    surface    reflected    waves. 

C.        WAVENUMBER    TECHNIQUE    APPLICATION 

Lauer  (1979)  recommends  a  receiver  that  is  a  single 
omnidirectional  hydrophone  and  a  source-generated  continuous 
wave  (Cil)  received  signal  for  source  depth  determination. 
He  further  s-ates  that  in  acuual  use  the  single  hydrophone 
could  be  replaced  by  an  array  ro  improve  the  signal-co-noise 
ratio  and   to   obtain    bearing  and    bearing  rate   information. 

An  example  of  a  possible  application  of  the  WT,  source 
depth  deTerminat ion,  is  presented  in  Fig.  3.  A  propagation 
model  could  be  run  for  a  series  of  incremental  source  depths 
•^0  yield  the  null  spacing  relationship,  as  was  shown  in  Fig. 
1.  Thai  curve  would  be  specific  to  -he  properties  of  the 
ocean  and  values  for  parameters  such  as  frequency  and 
receiver  depth.  A  received  signal  as  shown  on  the  righ-  side 
of    Fig.    3    could    then    be    analyzed.      lis   null    spacing   could    be 


28 


UT  APPLICATION:  SOURCE  DEPTH  DETERMINATION 
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Figure   3.      WT    Applica'iion    for   So  urea    Depth    Determination. 
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compared  to  the  model-producei  curve  to  infer   the  depth  of 
the  source. 

D.   TZCHNIQDE  LIMITATIONS 

The  success  of  the  WT  is  related  directly  to  the  accu- 
racy of  the  predicted  pressure  field  produced  by  the  propa- 
gation model.  There  are  several  possible  differences 
between  propagation  models  which  might  affect  the  pressure 
field.  Some  of  the  more  significant  include  approximations 
of  the  wave  equation,  bottom  effects,  initial  source  func- 
tions and  input  limitations.  These  will  be  compared  for  sev- 
eral models  to  determine  relevant  restrictions.  The  pressure 
field  will  also  respond  to  changes  in  the  ocean.  A  series 
of  different  scenarios  will  be  run  to  provide  an  initial 
esiimare  of  the  sensitivities  -o  some  of  these  factors. 

Any  additional  limitations  will  be  related  to  computer 
processing  time  and  storage  requirements.  These  factors  will 
not  be  addressed  directly  in  this  research  because  the  pro- 
cessing following  the  propagation  model  is  minimal.  There  is 
also  no  requirement  zc  store  the  complex  pressure  fields  or 
their  transforms  once  the  null  spacing  has  been  determined. 
Thus  the  WT  computer  'Jtilization  is  directly  proportional  "o 
the  time  and  storage  required  for  the  propagation  models. 
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III.   PROPAGATION  MODELS 

A.   INTRODUCTION 

Modeling  the  propagation  of  sound  in  the  ocean  is  com- 
plicated by  the  variability  of  the  ocean,  the  great  range  of 
frequencies  of  interest,  and  the  many  applications  of  sound 
propagation  (DiNapoli  and  Deaven port, 197  9)  .  It  is  not  sur- 
prising that  no  single  model  is  presently  capable  of 
addressing  all  of  these  variations;  there  are  many  semi-re- 
strictive methods.  In  general  all  models  consider  the  sound 
speed  to  be  a  function  of  the  spatial  coDrdinates  and  inde- 
pendent of  time  (DiNapoli  and  D eavenport, 1979) .  Only  the 
frequency,  water  column  georastry  and  oceanic  description 
will  affect  propagation. 

All  models,  regardless  of  application,  can  be  divided 
into  -WO  classes:  range  independent  or  range  dependent. 
Range  independent  models  assume  that  the  Dcean  is  cylindri- 
cally  symmetrical,  the  speed  of  sound  is  a  function  only  of 
depth  and  all  boundaries  are  planar  perpendicular  to  the 
depth  axis  (DiNapoli  and  Deavenpor-, 1979)  .  Range  dependent 
models  can  provide  better   apprDxima tions  of  real  conditions 
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by  allowing  the  sound  speed  to  be  a  function  of  two  or  three 
spatial  variables.  They  also  do  not  rsgiiire  planar  bounda- 
ries perpendicular  to  the  depth  axis,  thus  cylindrically 
syametric  bottom  topography  can  be  included  (DiNapoli  and 
Deavenport , 1979) .  A  range-dependent  modal,  the  split-step 
Fast  Fourier  Transform  PE  (SSFFT)  was  the  numerical  model 
used  in  this  research  to  generate  the  pressure  fields  for 
the  wavenumber  analyses.  This  PE  will  be  qualitatively  com- 
pared to  two  other  range-dependent  algorithms,  a  range-modi- 
fied Normal  Mode  (NM)  model  and  a  Finite  Difference  (FD) 
model.  It  will  also  be  compared  to  a  range-independent 
methcd,  the  Fast  Field  Program  (FFP)  moiel.  The  range-de- 
pendent capability  cf  the  SSFFT  was  not  used  in  this 
research  since  this  was  an  initial  investigation.  Its  com- 
parison with  both  range  dependent  and  independent  models 
may,  however,  indicate  possible  model  effects  on  the  WT 
which   could    be    studied  later. 

In  general  all  models  attempt  to  solve  the  acoustic  wave 
eguation.  This  equation  is  valid  for  sound  propagation  in 
fluid  media  with        specified  boundary  conditions 

(QricJc,!  975,p.  1 14)  .  It  becomes  the  rediced  elliptic  Helm- 
holtz  equation  in  spatial  coordinates  under  the  assumption 
cf    harmonic   time   dependence. 
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B.       PARABOLIC    EQUATION 
1  •      Introduction 

The  PE  algorithm  used  ia  this  rsssarch  was  divelop'Sd 
by      Brock    (1978).  The    basic      feature      of   this      PE   is      the 

replacement  of  the  recuced  elliptic  Helmholtz  equation  by  a 
parabolic  partial  differential  equation  that  utilizes  the 
Tappert- Hardin  split- step  Fourisr  algorithon  (SSFFT)  for  the 
numerical  integration  ( Brock, 19 78) .  The  original  computer 
code  as  described  by  Brock  has  been  modified  to  incorporate 
-he  source  and  a-tenuation  descriptions  ot  Tatro  (1977)  and 
to  include  an  ocean  bottom  by  Stieglitz  b^  al  (1979).  This 
is  the  version  of  the  ?E  that  Is  currently  available  at  NPS 
and  was  the  cede  (PEMODEL)  used  for  this  research.  (It  must 
be  emphasized  tha*  some  difficulty  was  encountered  in  trying 
to  ascertain  th?  correct  descr ip"^.ion  of  -he  NPS  version 
because  of  inadequate  and  erroneous  comment  card  documenta- 
tion in  the  source  cede  and  incomplete  manuscript  iescrip- 
tions.) 

This  section  will  discuss  t'he  SSFFT,  the  original 
Brock  algcrirhm,  modifications  to  -hat  cole,  and  subsequent 
•rvaluaticns  of  the  SSFFT.  Particular  a't-ention  will  be 
given   to    those   aspects  which   are    expected    tc    affect   the    WT. 
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2-      Parabolic   Approximation 

The  parabolic  approximation  to  ths  wave  equation  was 
first  introduced  by  Lsontovich  and  Fock  (1946)  as  a  solution 
for  the  propagation  of  electrotoa  gnetic  waves  along  the  sur- 
face of  the  earth.  Its  first  introduction  to  underwater 
acoustics  was  by  Hardin  and  Tappert  (1973)  who  used  the 
split-step  technique  in  conjunction  with  the  FFT  to  solve 
this  particular  approximation  to  the  wave  equation.  The 
result  was  the  SSFFT.  Tappert  (1977)  related  the  use  of  the 
parabolic  approximation  to  sound  channel  propagation  in  a 
waveguide  that  is  thin  vertically  and  elongated  horizontally 
to    the  range   of   the   first   convergence   zone   or   farther. 

Long-range  propagation,  with  which  the  SSFFT  is 
basically  concerned,  will  usually  be  for  low  frequencies 
because  volume  absorption  increases  strongly  with  frequency 
(Tappert , 1 977) .  The      maximum    elevation       (and      depression) 

angles  of  propagation  for  long  range  transmission  must  be 
small  zo  sa-.isfy  the  FE  (Tapper  t,  1977)  .  The  derivation  of 
the  SSFFT  will  be  discussed  later  in  a  sijiplified  version  as 
presented    by   Brock    (1S78). 

It  appears  that  the  SSFFT  is  intended  for  use  in 
situations    where      bottom    interaction    and      surface   scattering 
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would  be  negligible.  Thus  the  use  of  ttie  SSFFT  to  examine 
the  Lloyd's  Jlirror  effect  in  the  real  iforld  may  be  inappro- 
priate if  the  ocean  surface  is  sufficiently  rough;  addi- 
tional errors  may  occur  in  situations  where  bottom  effects 
are    significant. 

3-      Brock    Algorithm 
a.      Description 

The  SSFFT  is  valid  for  acoustic  pressure  in  a 
nedium  of  constant  density  with  a  monof requency  source  and 
cylindrical  symmetry  about  the  depth  axis  (Brock, 1978) .  The 
acoustic    pressure   can   be    written    as 


?(R,Z)     =    •J;(R,Z)    H    l(K    R)    . 

0  0'' 


where  the  zero-order  Hankel  function  of  the  first  kind 
relates  the  acoustic  pressure  to  an  outward  propagating 
cylindrical  wave  envelope  function  (Brock,  1 978)  .  This  is 
valid  because  of  the  assumption  that  at  lew  frequencies  all 
significant  energy  will  propagate  approximately  horizontally 
away  from  the  source  (Tappert ,  1  9  77)  .  The  asymptotic  form  of 
the    Hankel    function    is 

H    UK    R)     =    {2/ttK    R}''^    sxdCKK    R-tt/4)}       ,KR>>1, 
0  0  0  "  0  o  -' 
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if    the  receiver    is   many      wavelengths   from   the   source.        This 
leads   to     (Brock, 1978) 


The      additional     assumption     of      neglecting      the      far      field 
effects  , 


is      the   "parabolic"      approximation      for   radial      transmission 
(Brock, 1 978) .      The   result    is   the    Leontovich-Fock    (1946)    PE 


ii^    =    i(A+3)ij;  ^ 


where 


A    =    -    1 


2K       3z2 

0 


and 


3    =    Ko2    (n2    -    1) 
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which   is    the    form    used  by    Brocic     (1978). 
b.      Implementation 

Brock  (1978)  used  the  SSFFI  to  solve  the  PE 
because  it  has  several  significant  advantages  which  outweigh 
its  disadvantages.  Advantages  include  exponential  accuracy 
in  depth,  second-order  accuracy  in  range,  energy  conserva- 
tion, unconditional  stability  and  computational  efficiency. 
DisadvciQtages  are  a  uniform  mesh  and  periodic  boundary  con- 
ditions to  satisfy  the  FFT,  and  filtering  of  th^  sound  speed 
profile  by  smoothing  discontinuities  to  avoid  spurious  high 
angular- freguency  components  (Br ock, 1 973) .  Additionally  the 
algorithm  assumes  a  flat  pressur e-r=iease  surface,  a  vanish- 
ing field  at  the  maximum  depth  and  a  pseudo-radiation  condi- 
tion at  the  water-bottom  interface  by  smoothly  attenuating 
the    field     (Brock  ,1  978) . 

A   numerical    algorithm  of 


^(R+AR,Z)    =    e^^^^^"*"^^    -KR,Z)     , 


leads   to 


•MJl+AR,Z)    =    e^"^''^^      ?FT"^  {    e^^^^   FFT    (.,j;CR,Z))     } 
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or 


iARK    (a2-l)/2  iARK2/2K 

i(;(R+AR,2)     =    e  °  FFT"l{e  °FFT(  ii;(  R  ,Z)  )  } 


to  facilitate  the  computations  since  a  spatial  FFT  and  its 
inverse  are  required  to  transfDrm  between  i=pth  and  wavenum- 
ber  (Brock,  1 978)  .  This  procedure  is  implemented  by  two 
alternating  steps,  the  first  of  which  considers  propagation 
in  a  homogeneous  medium  to  account  for  diffraction  ard  the 
second  to  account  for  refraction  (Brock, 1 978) ,  In  his 
development  Brock  gives  an  alternate  calculation  of  the 
er.velope    function   stemming  from 

'MR-^AR,Z)     =    ei^^^/2    ^iARB    ^iARA/2       ^^^^^^^     ^ 


This  differs  slightly  from  the  previous  expression  contain- 
ing the  FFT  and  inverse  FFT.  Thus  it  is  :iot  currently  known 
which  relationship  is  used  by  Brock;  this  could  be  resolved 
later  by  a  detailed  examination  of  the  computer  code.  In 
any  event  as  McDaniel  (1975a)  has  shown,  either  approach 
would  be  of  sufficient  accuracy  with  appropriate  selection 
of    AH  and    AZ. 
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The  SSFFT  appears  to  be  gaite  satisfactory 
numerically,  but  it  has  to  simplify  the  actual  oceanic  vari- 
ability ro  allow  realistic  computation.  It  apparently  will 
handle  source  and  receiver  depths  and  horizontal  separation 
dis-ances  well  as  nored  by  its  degrees  of  accuracy  in  depth 
and  range.  The  smoothing  of  -he  sound  speed  gradients  and 
bottom  effec-^s    will   have    to  be   considered. 

The  intent  of  this  descriptiiDi  has  been  to  pro- 
vide a  base  from  which  comparisons  can  be  made  with  other 
numerical  solutions  and  modifications  to  this  algorithm. 
The  effect  of  the  source  model  on  the  WT  can  be  evaluated  by 
these  comparisons.  The  primary  elements  of  the  WT  are  the 
real  and  imaainary  parts  of  the  complex  pressure  amplitude. 
Thus  the  algorithmic  differences  with  resp  =  c-c  to  computation 
of  the  pressure  field  will  be  the  primary  basis  for  the  com- 
parison  of   effeccs   on   the    WT. 

McDaniel  (1975a)  compared  thr  basic  Tappert  and 
Hardin  SSFFT,  which  is  used  by  BrocK,  to  QDrmal  mode  theory. 
She  found  rhat  errors  arise  in  phase  and  group  speeds, 
although  one  mode  can  be  propagated  with  the  correct  ampli- 
tude and  group  ar.d  phase  spesis  .  McDanisl  describes  three 
sources    of   error:     (1)    approximating   the   wav?    equation    by   the 
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SSFFT  as  not^d  above,  (2)  limiting  the  range  steps  and  per- 
missible sound  speed  gradients,  and  (3)  truncating  the  field 
at  a  finite  depth.  The  result  of  the  SSFFT  will  be  a  shift 
in  the  lodal  interference  at  long  ranges  if  many  modes  are 
propagating  (McDaniel  ,1975a) .  The  error  due  to  range  step 
limitations  is  third  order  la  the  range  step  increment 
(McDaniel,  1975a) ,  thus  a  small  range  step  is  required  to 
minimize  this  error.  The  third  error  is  related  to  inclusion 
of  a  highly  absorbing  bottoi  which  causes  preferential 
attenuation  of  higher-order  modes;  the  dominant  contribution 
will   thus    result   from   the    lower    modes    (McDaniel, 1975a) . 

Brock  et  al  (1977)  have  described  a  procedure 
to  improve  this  condition.  It  maps  the  index  of  refraction 
and  depth  such  that  the  phase  speeds  and  turning  points  of 
the  ncrmal  modes  are  preserve!.  This  involves  construction 
of  a  "pseudoproblem"  such  that  the  parabolic  phase  speeds 
will  be  equal  to  the  elliptic  phase  speeds  of  the  corre- 
sponding modes  in  the  original  problem  (Brock  ^t  al  ,1977). 
This    is    giver,    as 


Cn,2)    -   (n,Z)    =    {(2n-l)^    ,    Zn^ 
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It  is  found  to  greatly  improve  the  SSFFT  agreement  with  the 
elliptic  solution.  This  corre::tion  technique  has  one  basic 
constraint,  which  is  that  isospeed  regions  are  not  allowed 
in  the  sound  speed  profile  because  of  poor  eigenf unction 
napping    (Brock   et   al,  1977)  . 

U.      Subsequent   Modifications 

a.  Introduction 

There  have  been  two  significant  modifications 
since  the  original  algorithm  was  published  by  Brocic  (1978). 
These  are  source  function  and  volume  attenuation  alterations 
(Tatro,1977)  and  the  incorporation  of  a  bottom  and  variable 
range  step  (Stieglitz  et  al  ,197  9).  The  Tatro  modification 
was  developed  in  1977  but  was  added  to  the  Brock  algorithm 
after  1978,  Discussion  of  these  two  aspects  will  complete 
the  description  of  the  PE  algorithm  installed  at  NFS  and 
utilized  in  this  study. 

b.  Tatro    Modifications 

The  Gaussian  source  originally  utilized  by  Brock 
(1978)  was  considered  to  be  inefficient  in  terms  of  program 
size  limitation  and  an  upper  frequency  limit  which  is  rela- 
tively low  (Tatro,  1977)  ,  a  new  source  function  was  devel- 
ooed.     This  was   intended  to   solve  the   problem  of   large 
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vertical  wavenumbers  from  a  Saiissian  sourcs  that  were  previ- 
ously  eliminated  with  range  (Tatro,1977|  and  led  to  an 
incorrect   depiction   of  the   field. 

This  source  is  basically  a  low  pass  filter  in 
*he  vertical  wavenumher  domain  which  is  centered  at  the 
source  depth  (Tatro, 1 977) .  The  initial  vertical  wavenumber 
field  will  thus  be  a  constant  u p  to  a  prescribed  value  and 
zero  for  higher  wavenumbers.  The  use  of  the  filter  is 
claimed  to  yield  improvement  d ver  the  Saussian  source  by 
minimizing  the  aliasing  which  might  occur  upon  transforming 
to  the  depth  domain  (Tatro, 1977|  .  The  benefits  derived  from 
the  use  of  this  new  source  are  reduced  computation  time  and 
storage  and  the  availability  of  higher  frequencies 
(Tatro, 1  977)  . 

He  also  modified  the  Brock  algorithm  by  intro- 
ducing volume  absorption  as  calculated  from  the  equation  of 
Thorp  (1967)  .  The  decision  can  be  made  whether  or  not  to 
include   volume   absorption. 

c.      Stieglitz    ejt    a2   Modifications 

Brock  (1978)  had  considered  the  ?E  most  useful 
for  long  range  propagation  a-  low  frequencies  and  small 
grazing  angles.    The   bortom   was   aodeled   as    fully    absorbing    so 
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that  the  field  would  vanish  at  the  maximum  depth  of  the 
traasfcrm.  This  depth  was  obtained  by  extending  the  water 
column  depth  by  one- fourth  to  obtain  the  transform  depth. 
The   index   of  refraction   in  this    region   was 


2  o  -{(Z-Z        )/b}2 

n^    =    n    2    +    ^^g  max    ' ""  ^ 


where  a  and  b  are  empirical  constants.  The  energy  in  the 
bottom  was  attenuated  such  that  any  additional  modes  which 
resulted    from      transform    truncation      were   minimized.  This 

empirical  expression  was  selec-ed  to  coincide  with  a  compa- 
rable normal  mode  solution.  This  rather  crude  model  of  the 
bottom  appeared  to  be  acceptable  for  situations  wherein  the 
bottom  had   a    relatively    small    =ffect. 

Stieglitz  et  al  (197  9)  considered  it  nscessary 
-0  attempt  to  develop  a  more  realistic  ocean  bot-om  since 
acoustic  propagation  was  being  utilized  in  more  varied  and 
extreme  situations  where  bottoa  interaction  was  important. 
They  developed  two  additional  options  for  handling  the  bot- 
tom: (1)  a  partially  absorbing  bottom  wi-h  specification  of 
reflective  less  versus  grazing  angle  of  the  equivalent  ray 
and     (2)         the    direcr    insertion      of   sediment    sound      speed   and 
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attenuation   profiles.  The   algorithm   derives      the  sediment 

profiles  from  the  grazing  angle  loss  values  in  the  first 
cpricn.  The  supplied  profiles  are  used  directly  in  the  sec- 
ond option  (Stieglitz  et  al  ,197  9).  A  new  complex  modified 
index  of   refraction 


M(R,Z)     =    {n2(R,Z)    -    1}    +    i2b(Z)/K 

0 


is    computed  if   the   bottom   attenuatiDn   is    to    be   calculated. 

One  additional  change,  the  calculation  of  the 
range  step,  has  also  been  addad  by  Stieglitz  et  al  (1979). 
I*:  consists  of  an  initial  search  of  the  vertical  wavenumber 
domain  for  the  maximuni  pressurs  component.  This  is  then  fol- 
lowed by  a  search  for  the  first  component  which  is  50  dB 
below     the      maximum      component.  The      wavenumber      of      this 

compcnen  t 


k      =    K      sine 

Z  O  50 


is    stated    to    reguire    a  range    step   of 


AR   =    X    /(I   -   cose      )  , 

0  5  0 


4U 


where 


X      =    2tt/K 

0  o 


The  range  step  was  fixed  in  the  current  research  for  the 
benefit  of  the  spatial  transform  in  the  wr,  thus  the  above 
description  was  not  utilized.  It  has  been  included  to  com- 
plete the  description  of  the  computer  cole  which  apparently 
comprises  the  NPS  model,  PEMODEL. 
5-       SSFFT    Evaluations 

The  SSFFT  has  been  investigated  dq  several  occasions 
to  examine  specific  items  which  might  be  modified  to  yield 
bet-er  agreemen-  with  actual  propagaxioa.  Most  of  these 
studies  deal  with  the  algorithmic  description  by  Tappert  and 
Hardin  (1973).  Examination  of  these  studies  will  note  fur- 
ther possible  limitations  on  accuracy  in  the  SSFFT  which  can 
be    extended   to   the    WT . 

Fitzgerald  (1975)  examined  the  accuracy  of  the  SSFFT 
as  a  function  of  range  and  frequency.  It  was  shown  that  the 
SSFFT  demonstrated  greater  accuracy  at  longer  ranges,  the 
lower     the    frequency.  The      SSFFT      was   considered      a      good 

approximation    in   the    examples   given    for      100    HZ    and    10    HZ    to 
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about  111  and  15000  km,  respectively  (Fitzgerald,  1975). 
The  current  study  will  examine  transmission  at  50  and  100  HZ 
at  ranges  to  100  nautical  miles.  These  ranges  will  exceed 
the    limits   of   Fitzgerald    in  some    cases. 

acDaniel  (197  5b)  investigated  splitting  the  total 
field  into  transmitted  and  reflscted  fields  in  deriving  the 
SSFFT.  It  was  necessary  to  neglect  the  reflected  field  for 
the  matrix  associated  with  the  Tappert-Hardin  SSFFT,  since 
the  two  fields  did  not  decouple  where  the  wavenuiiiber  was 
independent  of  range  (McDaniel,  1975b) .  Thus  an  exact 
expression  could  not  be  determined  for  the  transmitted  field 
without  the  reflected  field.  (The  reflected  field  is  neg- 
lected in  the  current  research  since  it  uses  the  SSFFT.) 
This  will  lead  to  phase  speed  and  group  speed  errors  as 
McDaniel  (1975a)  pointed  out  earlier,  as  discussed  in  Sec- 
tion III-B-U.  These  errors  were  minimized  by  Brock  et  al 
(1977)  by  a  mapping  technigue  for  both  the  index  of  refrac- 
tion and  the  depth.  Even  so,  the  implication  of  this  possi- 
ble error  en  the  WT  is  that  the  interference  pat-ems  may 
not  be  exact,  which  will  lead  to  inexact  or  biased  wavenum- 
ber  nulls.  McDaniel  (1975b)  proposed  another  matrix  which 
achieved   the    desired    decoupling    and   thus    allowed   expressions 


46 


for  the  transmitted  and  reflectad  fields.  This  alternative 
was  not  investigated  in  this  initial  research.  It  may  be  a 
way  to  resolve  non-exact  results  and  could  be  studied  later 
to  isolate  quantitatively  the  biases  related  to  various 
input   parameters   for    the    WT. 

Palmer  (1976)  related  possible  inaccuracies  in  the 
SSFFT  to  approximations  involving  relatively  small  terms  in 
the  exponential  approximations  for  the  operators  A  and  B  and 
also  to  the  use  of  a  single  reference  wavenumber.  These 
approximations  will  net  be  investigated  specifically  in  this 
research  but  they  may  be  related  to  potential  disagreement 
between  the  SSFFT  and  T:he  theoretical  Lloyd's  Mirror  effect. 
The  association  of  errors  due  to  a  singls  wavenumber  can  be 
recognized  from  the  wavenumber  spectra,  which  will  be  pre- 
sented later  in  section  IV-A.  They  will  show  not  only  an 
energy  concentration  at  the  reference  wavenumber  but  signif- 
icant energy  distributed,  as  well,  over  a  range  of  values. 
Palmer  (1976)  emphasized  that  the  physical  properties  of  the 
ocean  must  be  considered  before  terms  can  be  neglected,  such 
as    the   spreading   of   energy  over    several    wavenumbers. 

DeSantc  (1977)  studied  the  relationship  between  the 
acoustic    pressure   and  the   envelope   function    (called   velocity 
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potential  by  DeSanto),  which  is  directly  relevant  to  the 
current  research  where  pressure  is  calculated  from  the  mod- 
eled envelope  function.   DeSanto  arrives  at 


^'(R^Z)  =  A  R^'^^  /"?(t,Z)  Q(R,t,Z)  exp  (IK  /Zt)(R2  +  t2)  t^^^^^dt  ^ 


where  Q  satisfies 


QrR  ^  2{Cl-a)/R  +(iK^R)/t}Q^  +  Q^^  +  ^IK^Q^  +  2Q^l_{ln(p(t ,Z) ) } 

dz 


K^2{K2(t,Z)  -  K2(R,Z)}Q 


and 


K  (R,Z)  =  K(R,Z)  -  K.  (Z)^^ 


The  PE  for  the  envelope  function  is  obtained  if  this  func- 
tion is  evaluated  by  stationary  phase;  thus,  it  appears  zhat 
a  ncre  accurate  result  could  be  obtained  by  better  evalua- 
tion of  the  integral.  DeSanto  st  al  (1978)  developed  a  cor- 
rection for   the  SS??T,    known  as   the  Corrected   Parabolic 
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Approximation  (CPA) ,  that  handled  phase  inaccuracies  better 
than  the  SSFFT.  Further  study  might  addrass  the  response  of 
the  WT  to  the  CPA  with  respect  to  improved  phase  accuracy  as 
a  function  of  range. 

6«   Environmental  Sensitivity  of  the  PE 

This  final  section  in  the  analysis  of  the  PE  will 
consider  the  expected  environmental  sensitivity  of  the  PE, 
and,  by  extension,  "che  »T.  The  important  relationship  of  the 
PE  to  the  HT  is  the  input  of  the  complex  pressure  amplitude 
and  ultimately  the  resultant  wavenumber  distribution.  As 
discussed  earlier  in  Section  II- B,  the  wavenumber  could  be 
associated  to  the  source  angle  and,  by  extension,  to  the  ray 
type  or  family  (Lauer,  1979).  Consequently,  oceanic  parame- 
ters which  modify  the  transmission  of  carrain  rays  can  be 
related  zo   the  resulting  wavenumber  distribution. 

The  sound  speed  profile  will  favor  certain  ray  paths 
at  different  depths  in  the  water  column.  There  should  be  a 
corresponding  shift  in  ray  paths  and  hence  wavenumbers  as 
-he  source  and  receiver  depth  geometry  is  changed.  The 
direct  and  surface-reflected  rays  will  travel  at  different 
speeds  depending  upon  the  gradients  along  the  profile  and 
thus  lead  to  different   spatial  interference  patterns.   This 
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can  leal  to  wavenumber  shifts  and  changes  in  the  null  spac- 
ing. Profiles  used  in  the  current  study  will  include  real 
examples  and  isospeed  cases. 

Another  important  aspect  is  the  transmission  geom- 
etry. Generally  the  model  will  be  run  at  considerable  range 
and  in  relatively  deep  water,  which  will  be  consistent  with 
rhe  applicability  of  the  PE.  A  series  of  source  depths  with 
a  single  fixed  receiver  depth  will  be  run  to  generate  the 
wavenumber  null  spacing  values  that  will  formulate  a  deter- 
mination curve  for  each  case  of  geometrical  variation-  The 
distribution  of  these  resulting  curves  can  be  examined  to 
evaluate  the  effects  cf  transmission  gaomstry  on  the  WT. 

Two  frequencies  will  be  examined  in  a  similar  man- 
ner. This  would  be  one  of  the  most  important  parameters  to 
investigate  in  future  studies  since  it  would  have  a  direct 
bearing  on  operational  applications.  It  is  expected  that 
there  will  be  differences  in  the  resulting  wavenumber  spec- 
tra depending  on  the  input  frequency.  The  relationship 
between  frequency  and  transmission  geometry  will  also  be 
studied  since  these  are  the  easiest  paraaeters  to  vary  dur- 
ing an  experiment  and  will  lead  to  the  most  complete  analy- 
sis cf  model  variability. 
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The  surface  will  be  assumed  to  ba  planar  and  pres- 
sure-release so  that  there  will  be  complats  specular  reflec- 
tion. The  minimuiQ  source  depth  used  will  be  that  at  which  a 
definable  interference  pattern  3an  be  obssrved  for  the  spec- 
ified frequency.  The  slope  approaches  zero  near  the  surface 
on  the  source  depth  determination  curves.  This  result  may 
indicate  that  the  WT  in  its  current  form  is  probably  not 
applicable   for   sufficiently  shallow   sourcas. 

The  final  important  arsa  of  variability  will  be  the 
bottom  boundary  conditions.  There  are  several  options  now 
available  in  the  PEMODEL.  The  fully  absorbing  and  partially 
absorbing  bottoms  will  be  considered.  k  fully  absorbing 
bottom  will  be  used  in  most  cases  to  simplify  the  interfer- 
ence pattern  by  eliminating  any  rays  which  interact  with  the 
bottom.  Thus  the  basic  interference  pattern  may  be  evalu- 
ated. Various  bottoms  can  then  be  testsd  to  study  the 
effect  of  the  bottocp  on  the  *T.  Differing  bot-om  less 
curves  will  be  compared  in  a  fsw  cases  ro  demons-rate  quali- 
tatively that:  the  bottom  does  have  an  affrct.  The  third  bot- 
-om  cprion  of  a  realistic  sediment  sound  speed  could  also  be 
investigated  once  the  general  affect  of  the  bottom  has  been 
dr'ermined.      It      woulc  not    of far    any      additional   significant 
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results   at   zhts   point   beyond     those   realized    from  the    fully- 
and    partially-  absorbing    bottoms. 

All  input  parameters  are  related  and  thus  the  exami- 
nation of  any  one  particular  factor  must  be  considered  with 
respecr  to  the  full  range  of  variation  of  the  other  factors. 
The  number  of  combinations  of  parameters  would  be  very  large 
for  a  model  such  as  the  PEMODEL.  Thus  this  research  can 
only    begin      to   explore     general    relationships.  Additional 

research      is      necessary      to        delineate      specific     parameter 
impacts   on   the    WT. 

C.        FAST    FIELD    PROGRAM 

"'•      General    Descr  iption 

The  Fast  Field  Program  (FFP)  was  selected  as  the 
first  model  to  compare  to  the  PEMODEL  siacs  it  was  the  ori- 
gin of  the  WT  and  the  initial  testing  of  the  WT  by  Lauer 
(1979)  was  performed  with  the  FFP.  The  description  of  the 
FFP    which    follows   is    summarized    from    DiNapoli    (1971). 

The  basic  result  of  the  FFP  is  direct  numerical 
ir.-egra-icn  by  the  application  of  the  FFT  to  field  theory  to 
compute  propagation  predictions  in  a  minimum  computation 
time.      The    pressure   can    be   represented   in    the    form 

p(R,Z)     =    /F"/"    3(Z-,,Z„;K)    H    ^(KR)    K    dK     , 
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where  G,    the   Green's    function,   oust  satisfy 


d^G    +    {K^(Z)-k    2}G    =    -5(Z-Z    ) 
r  s 


and    the   associated   boundary  conditions.        The   field   equation 
can    then    be   written  as 

p(.Z,R^;     =vp    AKC2/TTi)''e       om/R        E         e         o 

m  mm  * 


where 


-      n  (  n      n         irNiz-'a         imR     AK 

=GCZ,Zo; K)Ke         o 
m  '    S      m        m 


In  this  form,  p  can  be  obtained  through  use  of  the  FFT  at 
each  range  increment.  Given  p  as  a  function  of  range,  the  WT 
could   be      applied   as    usual.  DiNapoli    compared      an    example 

from  the  FFP  with  the  normal  lode  prograi  of  Bartberger  and 
Aclcler  (1973)  and  found  excellent  agreement  between  the 
locations  of  the  peaks  and  tha  real  parts  of  the  eigenval- 
ues. The  only  significant  difference  is  that  the  FFP 
includes    more      higher   order   moiss     but   doss    not      include  the 
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first  four  modes.  DiNapoli  states  that  the  FFP  offers  a 
significant  reduction  in  computation  time  yet  still  provides 
reasonable  results  as  compared  to  normal  node  theory.  This 
is  accomplished  by  expressing  the  sound  velocity  profile  as 
one  or  more  exponential  functions  of  depth.  These  exponen- 
tial functions  allow  the  use  of  recurrence  relations  to 
quickly   calculate   the   input  to   the  FFT. 

The  FFP  offers  two  options  for  the  modeling  of  the 
bottom;  a  two-layered  fluid  bottom  with  specified  sound 
speeds  and  a  semi-infinite  bottom  of  one  sound  speed.  The 
propagation  effect  from  the  first  would  be  eventual  return 
of  rays  once  reflected  from  tha  second  bottom  interface.  The 
semi-infinite  bottom  would  not  allow  any  return  of  sediment 
rays. 

2.      Comparison   of   FFP    and    ?B 

The  FFP  is  a  range  inds pendent  mDdsl  and  the  PE  is 
range  dependent.  The  environment  was  considered  constant 
ever  the  range  in  the  current  research,  and  thus  the  two 
models  ran  be  compared.  The  FFP  uses  the  sound  speed  profile 
as  an  exponential  function  wnich  is  a  different  approach 
from  the  linear  segmentation  of  the  ?E.  The  number  of  expo- 
nential   functions     could    be      sufficient   -co      model    an      actual 
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profile  approximately  as  well  as  the  PE.  The  PE  bottom  mod- 
eling does  net  appear  to  offer  any  significant  advantage  in 
theory  over  the  FFP,  which  allows  for  a  semi-infinite  or 
layered  bottom.  It  is  likely  from  a  theoretical  standpoint 
that  the  PE  and  FFP  inputs  to  the  WT  could  be  similar  enough 
to  allow  comparison.  The  wavenamber  spectrum  shows  a  series 
of  peaks  which  result  from  the  nearby  singularities  of  the 
Greenes   function.  These   can      be   expressed      as   the      normal 

modes.  The  utility  of  the  WT  stems  from  the  spacing  of  the 
nulls   between    the  envelopes  which  contain    these    peaks. 

D.        NORMAL    MODE    TH20RI 

1 •      General    Descr  iption 

The  fundamental  importance  of  the  normal  mode  solu- 
tion, and  the  reason  that  is  it  used  to  check  the  validity 
of  other  models,  is  that  it  is  an  exact  solution  to  the  wave 
equation  (Kinsler  et  al  , 1982 ,p.  a30)  .  The  following  descrip- 
tion continues  from  Kinsler  et  il  (1982, pp. U30-U32) . 
The  solution    for    a   point    source    is    given   by 


{v2-32/c23t2)  p  =  1  $Xi-i^)     e^""^     , 


wner  6 


i      =  (0,Z<.)  . 
0        S 
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Analysis    yields 


p    =    e^*^"^    (-JTr)EZ^(Z<.)Z^(Z)H    Hk   R) 
n   n      b      n  o         n 


and    Zn  must   satisfy 


d^Z^    +    {a)2/c2(2)    -    K    2}z       =    0 
n  n        n 

dz2 


where  Kn  is  constant  and  appropriate  boundary  conditions. 
This  solution  is  for  trapped  modes  and  does  not  include 
those  which  are  evanescent.  Normal  mode  theory  can  be 
linked  to  the  WT  by  the  spatial  representation  of  the  wave- 
number 

K(Z)    =      u) 


C(Z) 


The  rays  in  'the  deep  sound  channel  which  have  the  same  local 
dirsction  of  propagation  correspond  to  the  normal  modes  by 
the  relationship 


cos  8  =  K-, 


K(Z) 
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which  extends  to  the  normal  mode  phase  speed 


C   =0) 
P    — 
K- 


where   1    designates  the  1th  mode. 

Normal  mode  theory  lends  itself  to  many  applications 
in  underwater  acoustics  and  it  may  be  formulated  for  both 
range  dependent  and  irdependent  situations.  A  dependent  form 
which  is  appropriate  for  comparison  with  the  SSFFT  was 
developed    by      Kanabis    (1975).  He   developed      a   normal-mode 

model  which  allows  large  changes  in  the  sDund  speed  profile, 
depth,  and  bcttom  composition  with  range.  The  range  domain 
is  segmented  into  regions  within  each  of  which  sound  speed 
and  bottom  composition  are  functions  of  drpth  but  not  range. 
The  normal  mode  solutions  for  each  region  are  matched  at  the 
interfaces. 

2.      Comparison   of    N!l    and    ?E 

It  would  be  expected  that  the  wr  would  respond  simi- 
larly to  ihese  two  models.  This  is  based  upon  this  example 
of  a  normal  mode  program  and  supported  by  the  previously 
discussed    evaluation    studies      involving    PS    and    NM.  The    NM 

solutions    are    generally   applicable   to   long    range   propagation 
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ar  lew  frequencies,  which  is  quite  similar  to  the  PE.  The 
NM  can  be  identified  closely  with  the  WT  in  a  physical 
sense,  because  of  the  interference  effects  in  the  NM  that 
result  from  several  multiple  reflections  (Officer,  1958, p. 
117).  Officer  further  states  that  the  nuaber  of  reflections 
will  increase  and  the  time  intervals  between  the  incidence 
of  successive  reflections  will  decrease  as  the  range 
increases.  This  might  be  comparable  to  depth  variations 
which  will  be  shown  for  the  WT.  The  NM  yields  exact  solu- 
tions for  modeling  with  planar  boundaries.  Thus  the  results 
of  a  WT  based  on  NM  would  be  limited  only  by  approximations 
which  were  intrinsic  to  the  WT.  A  comparison  of  PE-  and  NM- 
based  WT  data  should  aid  in  evaluating  the  impact  of  the 
particular    PE   algorithm. 

E.        FINITE    DIFFERENCE    METHODS 
"••      General    Descr  iptions 

One  final  group  of  methods  rfhich  has  evolved 
recently  are  the  finite  difference  methods  of  solution  for 
the  parabolic  equations  (Lee  and  Papadakis,  1979).  They 
seek  a  more  general  solution  which  will  be  appropriate  in 
shallow  water  or  where  bottom  interaction  will  be  important. 
Th9y      further    state      that   the      FD      and   ODE      methods   will      be 
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superior  to  SSFFT  because  an  artificial  bottom  is  not 
required  for  the  benefit  of  ths  FFT.  Explicit  and  implicit 
finite  difference  schemes  are  considered  as  well  as  an  ordi- 
nary differential  equation.  The  explicit  finite  difference 
(EFD)  scheme  is 


n+i     K  ^     n 
u      =  e   3   u 
m         —  m         ^ 

3R 


where   u      represents    a,  (B.n  ,Zm)  ,    Rq    is   the  nth   range   ooint,    and 
m 

2m    is  the   mth   depth      point.        The   implicit    finite   difference 
(IFD)    method    from   Lee   and    Botssas    (1982)    is 


-hK^ 


e         3R   u 
m 


u 


n 


m 


The    FD  scheme      is   consistent,      stable    and      convergent    in    the 
solution    of   the   ?£    (Lee    and  Papadakis,    1979). 
The   ODE   method  is 


d      u      =    a    u      +    bm 

—     .71  ni  ni       

dR  (AZ)2 


(u    ^,    -2u      +u      ^  )      , 
.11+1  m        m-l         f 


wh-2re    a   and    b   are   related    by 


u„    =    a(K    ,R,Z)u    +    b(K    ,R,Z)u^ 
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This   technique      is   stable,        CDasistent   aid      convergent   also 
which  supports      its   use.      One      additional   aspect   of     the   ODE 
approach    is  the   calculation  of  a    variable   range    step. 
2.      Comparison   of  FD    and   ODE    w^th    PE 

One  result  of  the  FD  and  ODE  nethods  is  greater 
flexibility  in  arbitrary  bottom  and  surface  boundary  condi- 
tions than  that  available  with  the  SSFFT.  Lee  and  Papadakis 
(1979)  state  that  while  the  FD  and  ODE  methods  are  bo^h 
superior  to  *he  SSFFT,  the  best  appears  to  be  the  ODE,  fol- 
lowed by  the  EFD ,  which  is  limited  by  step  size.  Thus  the  FD 
and  ODE  methods  appear  to  be  significant  improvements  over 
the  SSFFT.  This  result  should  extend  to  a  WT  based  upon  FD 
or    ODE   input. 

More  recently  Lee  et  al  (1981)  have  re-examined  the 
IFD  and  found  it  to  be  as  equally  efficient  as  the  ODE.  The 
general  advantage  of  the  IFD  over  the  SSFFT  is  that  the 
problem  is  solved  within  the  water  column  since  the  bottom 
effects  can  be  modeled  by  straight  line  segments.  This  is 
sifflilar  to  the  ODE,  tut  the  IFD  does  not  require  -.hs  amount 
of  storage  necessary  for  -^he  current  ODE  (Lee  et  al,  1981). 
McDaniel  and  Lee  (1982)  have  extended  the  IFD  to  the  treat- 
ment  of      ver"ical    density    discontinuities.       They      found    that 
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while  without  this  more  realistic  treatment  of  the  water- 
bot-om  interface  the  IFD  approximates  the  SSFFT  and  differs 
significantly  from  normal  mode,  the  IFD  corresponds 
extremely  well  with  the  normal  mode  calculations  with  the 
interface.  In  conclasion,  the  numerical  techniques  related 
to  FD  and  ODE  appear  to  be  viable  alternatives  to  the  SSFFT. 
They  will  permit  a  more  accurate  representation  of  underwa- 
ter acoustic  transmission,  which  will  be  directly  relatable 
to    the  WT. 
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IV.   ANALYSIS  OP  WT 

A.   GENERAL  DZSCBIPTION 

Twenty  scenarios  were  sslscted  to  examine  the  general 
aspects  of  variations  in  several  oceanic  and  geometric 
paraaeters.  A  set  of  input  parameters  and  a  sound  speed 
profile  were  constant  in  each  scenario  while  the  source 
depth  was  varied.  Five  x.o  ten  source  depths  were  used  for 
each  situation.  The  vertical  and  horizontal  axes  of  the 
wavenumber  spectra  correspond  to  the  spectral  intensity  and 
■^he  scaled  wavenumber  beta,  respectively.  The  intensity 
axis  was  normalized  to  unity  for  ease  in  plotting.  The  null 
spacing  was  determined  for  each  spectrum.  These  null  spac- 
ing values  were  plotted  versus  source  depths  to  yield  the 
determination  curve  for  each  scenario. 

Table  I  lists  the  parametric  specifications  for  each 
scenario.  Basically  they  car.  be  divided  into  isospeed  or 
sound  speed  profile  cases.  An  isospeed  profile  of  4890  ft/s 
from  the  surface  to  the  bottom  was  sei  =  c'-ed  *■  o  evaluate  the 
transmission  geometry  and  the  general  consistency  of  experi- 
mental results,   assuming  a  fully-absorbing  bottom.    The  WT 
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is  based  on  rhe  arrivals  of  the  direct  and  surface-reflected 
rays;  thus,  the  isospeed  cases  will  be  the  easiest  to  com- 
pare with  theory.  There  will  be  no  sound  speed  differences 
as  a  function  of  depth  which  could  modify  the  arrival  times. 
Additionally,  several  realistic  ocean  profiles  were  studied 
to  allow  a  qualitative  comparison  with  theory  and  explore 
spa-ial  and  temporal  variability.  It  is  recognized  that  the 
rheoretical  Lloyd's  I!irror  sffecr  applias  only  to  isospeed 
profiles.  Results  from  real  profiles  were  compared  with  this 
theory  zo  see  if  they  could  approximate  it.  The  spatial  and 
tempcral  scales  of  variability  used  in  this  research  were 
approxiaately    200    miles   and  20  days. 

More  analyses  were  performed  with  "rhe  isospeed  cases 
since  this  was  an  initial  investigation.  It  still  remains 
to  consider  other  simple  cases,  such  as  Linear  negative  and 
positive  gradients,  surface  ducts  and  deep  sound  channels. 
Also  realistic  ocean  studies  that  would  involve  greater  tem- 
poral and  spatial  variations  wDuld  be  beneficial.  Time  was 
limited  during  this  research,  thus  efforts  were  divided 
between   model      validation    and      oceanic   sensitivity.  It    is 

believed,  however,  that  the  comparisons  which  will  be 
presented  are  sufficient  for  a  preliminary  investigation 
into   the    sensitivity    cf    the   WT. 
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TABLE    I 
Scenario  Specifications 


ISCSPEED    PROFILE 


Run 
Set 

Pcv  Depth 
RD:ft 

Range 
R  :nm 

Freg 

F:az 

Water  Depth 
Z:ft 

Number  of 
SD  Runs 

10 

300 

50 

100 

10003 

7 

11 

300 

50 

50 

10000 

6 

12 

300 

100 

50 

10000 

6 

13 

300 

50 

100 

5000 

6 

lU 

300 

100 

100 

5003 

5 

15 

300 

50 

50 

5000 

6 

16 

300 

100 

50 

10003 

6 

17 

30  0 

100 

100 

10003 

5 

18 

800 

50 

50 

10003 

5 

19 

10000 

50 

50 

10003 

5 

Rl 

300 

25 

100 

10000 

5 

22 

300 

50 

100 

10003 

5 

R3 

300 

75 

130 

10003 

5 

REALISTIC    PROFILE 


Run 
Set 

Rev  De 
RD:f-^ 

pth 

Range 
R  :nm 

Fr-g 
F:Hz 

Water  Depth 
Z:ft 

Src 
SD(i 

Depth 
lumber) 

AO 

300 

100 

100 

9000 

8 

A1 

1000 

100 

100 

9003 

5 

A2 

9000 

100 

100 

9003 

5 

A3 

300 

100 

100 

9003 

5 

50 

30  0 

100 

130 

12003 

8 

CO 

300 

100 

100 

9003 

8 

DO 

300 

100 

103 

12003 

8 

ea 


Th€  complex  pressure  input  data  wers  not  modified  so 
that  tha  entire  waven  umber  spectrum  would  be  observed  at  the 
receiver.  These  effects  would  need  to  be  evaluated  if  the 
WT  was  to  be  compared  to  received  signals.  The  surface  and 
bottom  were  assumed  flat.  Boundary  variations  would  have  to 
be  considered  for  WT  use  under  realistic  ocean  conditions. 
A  nominal  beam  width  of  thirty  degrees  was  used  and  in 
almost  all  cases  the  range  step  was  0.1  mile. 

Frequencies  of  50  and  100  Rz,  water  column  depths  of 
5000  and  10,000  feet  and  ranges  of  50  and  100  nautical  miles 
were  the  basic  comparisons  fDr  most  isospeed  scenarios.  A 
receiver  dep-^h  of  300  feet  and  a  fully  absorbing  bottom  were 
utilized,  unless  otherwise  noted.  Additional  comparisons 
were  performed,  such  as  range  variations  to  25  and  75  miles 
and  receiver  depth  variations  of  800  and  10,000  feet.  A 
range  of  100  nautical  miles  and  a  frequency  of  100  Hz  were 
used  for  the  aciual  cases  unless  otherwise  noted.  Th=  bottom 
depths  were  9000  feet  for  profiles  A  and  C  and  12,000  feet 
for  profiles  B  and  D. 

The  realistic  profiles  were  obtained  from  the  Acoustic 
Storm  Transfer  and  Response  Experiment  (AST3EX)  which  was 
conducted  by  NPS  in  the   norrheast  Pacific  Ocean  in  November 
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and  December,  1980  (D unlap,1982)  .  The  ASIREX  study  was  con- 
ducted to  evaluate  the  effect  of  oceanic  response  to  atmos- 
pheric storms  or.  acoustic  propagation.  The  profiles  were 
selected  for  this  research  zo  represent  only  spatial  and 
temporal  variability  without  respect  to  any  experiment  oper- 
ations. The  profiles  were  obtained  from  expendable  bathy- 
thermographs (AXBT)  which  were  dropped  by  U.  S.  Navy  P-3 
aircraft.  Fig.  U  shows  the  relative  locations  and  times  for 
these  profiles.  The  digitized  profiles  were  analyzed  by  the 
University  of  Hawaii  to  remove  spikes  and  false  values. 
These  profiles  were  combined  with  climat ological  salinity 
profiles  at  the  Fleet  Numerical  Oceanography  Center  in  Mon- 
terey. Table  II  gives  the  descriptive  data  for  each 
profile.  The  bottom  loss  curves  which  were  used  are  similar 
to  those  found  in  Orick  (1979)  and  were  arbitrarily  devel- 
oped to  indicate  higher  loss  for  profiles  B  and  D.  This  is 
related  to  bottom  composition;  bottom  types  for  profiles  A/C 
ar.d  3/D  are  calcareous  sand  and  clay,  respectively 
(NA70CFANO,1 978)  .    Thus  higher    loss   is    expec'ted   for    B/D. 

A  general  feature  of  the  spectra  for  almost  all  rest 
cases  is  a  distinct  (J-shape  with  peak  values  at  the  minimum 
and    maximum    wavenumbers.       Figs.       5   through    10    show   this    phe- 
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TABLE    II 
Location   Information   for    Realistic   Profiles 

Profile        Lat  Long       Date        Time(Z)       MLD  (m) 

A  U1-2aN  127-23W  16N0V  1928:0a  32 

B  U5-a7N  133-05W  15N0V  2121:25  61 

C  U1-23N  128-22W  6  DEC  18U3:37  58 

D  45-46N  133-09W  5  DEC  2033:45  79 

r.oaenon   for      isospeed  cases  and    Figs.         11    through    16      for   a 
realistic   sound   speed   profile.  The   abbreviations   on    these 

figures  are  source  depth  (SD)  ,  receiver  depth  (RD)  ,  and 
water  column  dep-h  (Z)  in  feet;  range  (R)  in  nautical  miles; 
and  frequency  (F)  in  Hz.  This  key  will  also  apply  to  the 
source  depth  determination  curves  which  will  be  presented 
later.  The  righ-^  and  left  intensity  maxiia  appear  lO  corre- 
spond to  beam  eleveation  angles  of  0  and  30  deg,  respec- 
tively, which  could  be  considered  to  approximate  the  direct 
and  surface-reflected  waves.  k  30  deg  b=am  angle  was  used; 
thus,  there  may  be  a  relationship  with  the  left  maximum. 
This  beam  angle  is  also  near  the  maximum  angle  which  Tappert 
(1977)  considered  to  be  appropriate.  Then  the  l^ft  maximum 
might  be  a  result  of  algorithmic  inaccuracy.  Different  beam 
angles  could      be  tested    to     determine    if    this      would    explain 
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th€    unexpected      0-shape   curvs.  This   apparently     anomalous 

spectral  pattern  might  be  explained  by  a  possible  transform 
assumption  of  symmetry  about  zaro  and  ths  inclusion  of  neg- 
ative Transform  values.  It  is  also  possible  that  the  left 
wavenumber  peak  may  be  the  result  of  aliasing  of  the  high 
wavenumber  values.  DiNapcli  (1971)  suggested  that  aliasing 
could  occur  ?nd  thus  he  utilized  only  ana-half  of  the  range 
field.  This  has  not  been  tested  in  the  current  research; 
however,  examination  of  different  range  segments  might 
resolve  xhis  apparent  problem.  A  remots  possibility  is  that 
the  direct  and  surface-reflected  rays  eaoh  exhibit  a  Bessel 
function  display.  This  might  be  depicted  by  a  peak  at  the 
reference  wavenumber  for  each  ray  with  a  set  of  interfering 
envelopes  be-cween  the  peaks.  A  marhemaci-al  explanation  for 
this  possibility  has  not  yet  been  explored.  No  spikes  ar=^ 
seen  in  very  shallow  vater  as  iepictsd  in  Fig.  M,  but  also 
no  dis-inct  null  pattern  is  observed,  thus  near  surface 
source   depths    have   be<=n    neglected. 

These  figures  were  also  included  to  show  a  series  of 
source  depth  runs  and  the  resulting  decrease  in  beta  spacing 
as    a    funcrion    of   source    depth.  The    aull    -pacing    was    manu- 

ally   measured    for   each  run   as    the   distance    between   two    adja- 
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Figure  5.      Iscsosed   Profile   Wave  number  Sorctrum    for    SD    of 
200      ft. 
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RUN  145 
SD:500FT  RD:300FT  F:100HZ  Z:5000FT  R:10 
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Figure   7.      Iscspeed    Profile   Wavanumber  Spactrum    for    SD    of 
500    ft. 
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Figure  9.   Isc  =  peed  Profile  Wave  numb 
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Figure   10. 


IscsDeed    Profile    Wavenumber   Spectrum    for   SD   of 
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RUN  R03 
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RUN  R04 
SD:300FT  RD:300FT  F:100HZ  Z:9000FT  R:100NM 


0.050   a.Oii       0.040       0.035       0.030       0.025       0.020       0.035       0.010       0.005       0.000 

SCALED  WflVENUMBER  BETA    11 /FT J 


Figure    12.      Realistic   Profile   Wavenumber   Soectrum    for    SD   of 
300    ft. 
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Realistic   Profile    Wavenmnber   Spectrum    for    SD   of 
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Realistic   Profile    Wavenumber   Soectrum    for    SD   cf 
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RUN  fl08 
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Figure   16. 


Realistic   Profile   Wavenumber   Spectrum    for    SD   of 
3C00   ft. 
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RUN  141 
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Fiaure    17.      Isosoeed    Profile    Wavenumber   Sosctrum    for    SD    of 

^  Qn    -Si 


82 


csnt  intensi-ty  minima,  such  as  shown  in  Figs.  5  through  16. 
The  nulls  selected  for  measurement  were  chose  in  the  middle 
of  the  spectrum,  which  led  to  a  consistent  technique  for 
comparison  and  avoided  the  rapidly  decreasing  null  spacing 
near   the    left   intensity    maxima.  The   distance   was   measured 

to  the  nearest  1/64  of  an  inch  and  equated  to  the  respective 
null   spacing.  Thus,      there      is   some      iiheren-    measurement 

error.  The  exac-  null  locations  were  difficult  to  determine 
in  the  shallow  isospeed  cases  because  of  gen-^  le  spectrum 
curvature  and  uncertainty  in  null  location.  The  actual 
profiles  were  very  complex  and  again  the  nulls  were  diffi- 
cult to  locate.  One  value  for  each  spectrum  plot  was  usad 
*o  generate  a  point  on  the  determination  curve,  such  as  Fig. 
19  for  a  series  of  runs.  The  variation  of  the  data  poin-s 
from  the  theoretical  Lloyd's  Mirror  curve  will  be  the  normal 
distance  between  "he  point  and  the  curve.  Measurement  inac- 
curacy will  decrease  as  source  depth  incraases.  This  will 
be  significant  because  at  shallow  depth  there  is  some  varia- 
tion noted  from  sound  speed  differences.  The  nulls  are 
well-defined  at  areat  depths,  but  there  is  lit-le  varia-^ion 
already,  thus  the  measurement  errors  were  net  considered 
significant.    It    is   possible   that    a   numerical    algorithm    which 
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could  quanti-atively  identify  the  null  locations  could  mini- 
mize  this   error. 

This  expectation  is  inferred  from  ths  wavenumber  plots 
by  Lauer  (1979)  and  his  examination  of  range  determination. 
He  found  tha-^  range  determination  resolution  increased  as 
the  distance  between  source  and  receiver  decreased.  It  is 
possible  that  a  similar  relationship  might  exist  for  depth 
dererminatior.  Thus  greater  resolution  aay  be  realized  at 
shallow  source  depth.  Even  at  greater  source  depths  the  res- 
olution may  be  acceptable  for  certain  applications  where  an 
approximate   depth   is    acceptable. 

The  prediction  curves  which  result  from  the  individual 
WT  spectra  will  be  discussed  in  following  sections.  Table 
III  is  a  listing  of  the  prediction  curve  comparisons  which 
were  conducted.  In  the  figures  only  on?  symbol  may  appear 
at  a  source  dep-h  if  both  curves  had  the  same  null  spacing 
value. 
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SOURCE  DEPTH  DETERMINATION  CURVE 
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Figure    18,      ::;xampl9   of  Source   Depth   Detarmination   Curvi 
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TABLE    III 
Comparison  Spscifications 

ISOSPEED    PROFILE 

Figure      Run   Sets  Constant  Parameters  7arying   Parameter 

19  11,15  R: 50,F:50, RD:  300  Z:10000,5000    ft 

20  10,13  R: 50,F: 100,RD :300  Z:10000,5000    ft 

21  12,16  R:  100,F:50,RD  :300  Z:10000,5000    ft 

22  17, lU  R:  100,F:100,  RD:300  2:10000,5000    ft 

23  15,16  F: 50,Z:5000,RD:300  R:50,100    nm 
2U             11,12  F: 50, Z: 10000, RD:300  R:50,100    nm 

25  R1,R2,R3  F :  100 , Z: 10003 ,  rD: 300  R:25,50,75    nm 

26  12,17  R:100,Z:10000,RD:300  F:  100,50     hZ 

27  10,11  R: 50, Z: 10000,  RD:300  F:100,50    hz 

28  11.18,19  R: 50, F:50, Z:  1 0000  RD: 300 , 800, 10000    ft 

REALISTIC    PROFILE 

Figure      Run    Sets  Constant   Parameters  Varying    Parameter 


29  AO,BO  R:  100,F:100, RD:300  Space:200    nm 

30  CO, DO  R:  100,F: 100, RD: 300  Space:200    nn 

31  AO,CO  R:  100, F:100,RD:300  rim3:20    days 

32  EC, DO  R:  100,F:100,rD:300  rime:20    days 

33  A0,A3  F:  100,  16N0V,  RD:300  R:100,50    nm 

3U  A0,A1,A2  F:  100,  16N37,   R:10Q  RD :  300  ,  1000, 9000    f- 

COMPOSITE 

Figure  Run    Sets  Constant  Parameters  Varying    Parameter 


35       (All   Isospeed  and   Ac-ua 1   Profile    Run    Sets) 
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B.        ISOSPEED    CASES 

1*      Water   Column    Ce^th  Varia  tion 

Fig.  19  is  the  comparison  for  a  frequency  of  50  Hz 
and  a  range  of  50  nm  at  water  column  iepths  of  5000  and 
10,000  feet.  It  was  not  expacted  that  there  would  be  any 
variation  if  the  water  column  depth  was  modified,  since  all 
bottom      interacting   rays      should      be      fully   absorbed.  The 

curves  showed  excellent  agreement  at  source  depths  of  800 
feet  and  greater,  but  there  was  some  variation  noted  for  500 
feet  and  shallower.  Some  variation  may  be  related  to  meas- 
uremenr      error,      as      noted   abov9  ,  Maximum    variation      from 

*heory  is  observed  for  both  iept-hs  in  the  middle  depth 
range,  approximately  500  feet,  which  is  also  the  region  of 
greatest    curvature    in   the    sourca    depth   determination   curve. 

Fig.  20  is  the  same  except  the  frequency  has  been 
increased  to  100  Hz.  There  is  better  agreement  from  200  to 
3000  feet;  however,  the  significance  of  the  frequency  varia- 
tion has  not  been  evaluated.  Fig.  21  is  for  a  frequency  of 
50  Hz  and  a  range  of  100  nm,  and  Fig.  22  is  for  100  Hz  and 
100  nm.  Some  variations  are  noted;  however,  they  may  not  be 
significant,      which      would  agree      with   ths      supposition   that 
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Figure   19. 


Iso^DG^d    Profile   SD    Detsrmination   Curve:    Z 
Vana-icn   -    1. 
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water  column  depth   differences  will  not  be   significant  for 
an  isospeed  profile. 

2.   Range  Variation 

Fig.  23  is  the  comparison  for  a  frequency  of  50  Hz 
and  a  water  column  depth  of  5000  feet  for  50  and  100  nm. 
There  is  essentially  no  difference  between  the  curves.  The 
experimental  curves  agree  well  from  800  to  3000  feet,  but 
once  again  the  shallow  variations  may  result  from  inaccura- 
cies in  evaluating  beta  from  the  spectrum  plots.  Fig.  2U  is 
a  similar  comparison,  but  the  water  column  depth  has  been 
increased  to  10,000  feet.  Fig.  25  increases  the  frequency 
to  100  Hz  for  ranges  of  25,  5  0  and  75  miles.  Again  all 
curves  show  very  good  agreement  with  theory,  which  indicates 
tha-  range  variations  within  103  nm  of  the  source  will  prob- 
ably be  insignificant.  Ranges  beyond  100  nm  were  not  tested 
since  Fitzgerald  (1975)  suggested  that  the  SSFFT  was  not  a 
good  approximation  for  100  HZ  beyond  111  ki.  This  could  be 
tested  later. 

3«   Frequer.cy  Variation 

Fig.  26  is  the  comparison  for  a  range  of  50  nm  and  a 
water  column  depth  of  5000  feet  for  50  and  100  Hz.  Both 
curves  agree   well  with  theory  below   800  feet,   but   at  500 


89 


SOURCE  DEPTH  DETERMINRTION  CURVE 
ISOSPEED  WRTER  DEPTH  VRRIRTIQN 
R:50NM  F2IOOHZ  RD:300FT 


LJ 
O 


LEGEND 
•  -  LLOYDS  MIRROR 
■  -  10:  Z-IOOOQFT 
a  -  13:  Z-5000FT 


0.0000  0.0025    O.OOSO    Q.OOTS    O.OIOQ    0.012S 

BCTR  NULL  SPflClNC   (l/rT) 


T" 


T" 


n 


o.oiso      i).oi;rs      0.0200 


Figure  20. 


Iso^oeed    Profile 
Variation   -    2. 


SD    Determination   Curve: 


90 
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faet  the  higher  frequency  appears  to  agree  better.  Fig.  27 
extends  the  water  column  depth  to  10,000  feet  and  there  is 
apparently  better  agreement  with  the  theoretical  curve. 
Overall  there  appears  to  be  better  agreeient  at  100  Hz  than 
at  50  Hz,  but  this  result  is  not  certain  with  the  limited 
number  of   cases   run. 

**•      Receiver   De£th  Yariat  ion 

Fig.  28  compares  receiver  depths  of  300,  800  and 
10,000  feet  for  a  range  of  50  nm,  water  column  depth  of 
10,000  feet,  and  frequency  of  50  Hz.  Again  all  curves  agree 
reasonably  well.  There  is  perhaps  slightly  better  agreement 
for  the  shallow  receiver,  but  this  supposition  cannot  be 
verified    from   these   limited  data, 

5.       Analysis   of    Isqspeed   Cases 

There  do  not  appear  to  be  any  significant  differ- 
ences for  warer  column  depth,  range,  fra^uency  or  receiver 
depth  over  the  limits  of  variations  considered.  This  was 
expected  since  the  theory  is  supposed  -o  be  only  a  function 
of  source  depth  for  ar  isospeed  profile.  Thus  -he  WT  would 
be  expected  to  be  equally  effective  for  various  transmis- 
sion geometry  patterns  and  frequencies,  if  the  sound  speed 
profile    is    constant    and    the   bottom   fully    absorptive.         There 
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may  be  slight  effects  for  receiver  depth  and  frequency,  but 
these   have  not   been   substantiated. 

C.        REALISTIC    PROFILE    CASES 
"'•      Spatial   Variation 

Fig.  29  is  th?  comparison  of  spatial  variability  for 
the  curves  which  resulted  from  the  spectra  of  profiles  A  and 
B.  (The      sound  speed      profiles   corresponding      zo  the      pre- 

diction curves  will  be  included  for  each  comparison.)  These 
were  both  sampled  on  November  16,  1980,  and  represent  a  spa- 
tial separation  of  abcut  200  nm.  Both  profiles  result  in  a 
slight  bias  below  the  theoretical  Lloyd's  Mirror  curve. 
There  is  more  variation  along  each  curve,  which  is  to  be 
expected  because  the  profiles  are  not  isospeed.  Also,  par- 
tially absorbing  bottoms  are  now  used  as  compared  to  the 
fully  absorbing  bottoms  used  with  the  isDspeed  cases.  The 
profile  and  bottom  effects  have  not  yet  bBer.  separated  to 
de-ermine  which  contributes  lore  to  -he  departure  from 
theory.  This  would  require  an  extensive  sBt  of  intercompari- 
sons   between    these   two  paramet  =  rs.  (The    present   intent    is 

•o  show  rhat  ocean  variability  is  important.)  It  is  diffi- 
cult to  dif f eren-ia-'-e  between  -he  -wo  curves,  but  it  appears 
•here   is      a      slight    trend      in      3      -.3 ward   smaller      null 


-ha- 
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spacings.  Sound  speed  profila  B  is  slowar  than  A  over  the 
entire  depth.  A  slower  sound  speed  would  appear  to  indicate 
that  a  ray  of  lesser  angle  would  be  required  to  provide  the 
second  arrival  for  the  interference  pattern.  That  would  cor- 
respond to  a  greater  vertical  wavenumber  and  decreased  beta 
null   spacing. 

Fig.  30  compares  profiles  C  and  Q,  which  were  sam- 
pled at  the  same  locations  as  profiles  A  and  B,  respec- 
tively, but  on  December  6,  1989.  A  similar  pattern  is  noted 
with  D  spacing  less  than  C,  which  may  again  correspond  to 
slower  sound  speeds  fcr  profile  D. 
2.      Temporal   Variation 

Fig.  31  compares  profiles  A  and  C  which  represent 
the  same  location,  tut  a  separation  of  20  days.  There 
appears  to  be  a  curve  crossing  at  about  800  feet,  such  that 
spacing  is  less  at  A  for  shallow  dep-chs  and  greater  for  A 
below  80  0  feet.  The  A  sound  sppsd  profile  is  slower  than  the 
C  profile  down  to  800  feet  and  the  pattern  is  reversed  from 
800  to  2C00  feet.  Slower  speeds  appear  to  be  related  to 
shor-er  null  spacings,  which  appears  coasis-cent  with  this 
example.  The  variations  between  the  curves,  however,  are 
quite    small. 
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Fig.  32  presents  a  similar  comparison  for  profiles  B 
and  D.  The  speed  profiles  show  almost  hd  differences,  except 
-hat      D      is   slightly      slower.  This      result   appears      as      a 

slightly    sharper   deep   sound  channel   which    will   mean   slightly 
slower  velocity      for    the      direct    path.         That   effect      may    be 
reflected    as   a   slightly   greater    null    spacing   as   shown. 
3.      Range    Variation 

Eesults  for  profile  A  wsr«=^  compared  for  ranges  of  50 
and    100      nm   in      Fig.       33.  No    significant      variations    were 

noted,  which  is  similar  to  the  results  for  -^he  isospeed 
cases. 

**•      Receiver   De£th  Variation 

Results  for  profile  A  were  also  compared  for 
receiver  depths  of  300,  1000  and  9000  feet  in  Fig.  34. 
There  may  be  slight ly  better  agreement  for  the  shallow 
receiver,    but    this   is  not    certain. 

D.       COMPOSITE    COMPARISON 

Fig.         35    is      a      composite   of     all      isospeed   and      actual 
profile      data    sets.  Isospeei      data      show    slightly      better 

agreement  wi-h  theory;  however  realistic  profiles  are  a  rea- 
sonable approximation  to  the  theoretical  curve.  Realistic 
profile      variations   from      the    isospeed      curve      appear    to      be 
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related  tc  the  sound  speed  profile  and  different  bottom 
boundary  conditions.  The  general  relationships  of  the  data 
indicate  that  below  about  800  fset  no  sigaificant  variations 
are  noted.  Specific  curves  will  apparently  provide  better 
predictions  at  shallower  depths.  It  may  be  sufficient  for 
certain  applications  zo  provide  climatological  or  theoreti- 
cal prediction  curves.  These  curves  coald  be  enhanced  or 
modified  by  updated  curves  that  would  reflect  the  current 
sound   speed   profile. 
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V.   CONCia SIGNS 

The  Wavenumber  Technique  ( WT) ,  as  proposed  by  Lauer 
(1979)  has  been  investigated  for  source  depth  determination. 
Several  idealized  and  actual  scanarios  of  oceanic  and  acous- 
tic parametric  variations  vera  tested.  There  is  excellent 
agreement  between  the  ideal  test  cases  and  the  theoretical 
Lloyd's  Mirror  effect.  The  raalistic  rases  showed  signifi- 
cant similarity  as  well,  but  with  some  bias  toward  smaller 
wavenumber  spacing.  This  is  qualitatively  attributed  to 
sound  speed  profile  variations  in  the  watsr  and  sediment. 

There  appears  to  be  a  response  of  tha  WT  to  variations 
in  input,  which  is  beneficial  in  allowing  evaluation  of  the 
relative  importance  of  individual  parameters.  Changes  in 
range  and  water  column  depth  do  not  indicate  any  significant 
variation.  Some  differences  are  possibly  suggested  for 
receiver  depth  and  frequency,  but  the  results  have  not  been 
fully  analyzed.  The  »?T  appears  to  respond  well  to  changes 
in  the  environmental  parameters,  su::n  as  sound  speed 
profile.  This  is  seen  from  the  average  bias  cf  the  actual 
cases  from  the   isospeed  cases.   The  bias  does   appear  to  be 
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rslativsly  consistent  and  shows  the  importance  of  oceanic 
variations  in  both  space  and  in  time.  Thus,  it  appears  that, 
if  information  is  available  for  the  bottom  structure,  and  if 
a  realistic  sound  speed  profile  can  be  generated,  then  a 
very  specific  family  cf  curves  could  be  generated  for  vari- 
ous frequency  and  receiver  depth  combinations.  This  simplic- 
ity indicates  possible  operational  applicability  and 
flexibility. 

These  are  the  first  known  source  depth  determination 
curves  for  the  WT  and  thus  there  is  no  couparison  available. 
It  may  be  that  there  are  some  biases  in  these  curves  which 
would  ba  discovered  *ihen  comparative  runs  are  made  using 
other  inputs  or  propagation  models.  Comparison  with  observed 
signals  might  also  reveal  biases  that  are  related  to  the  '^T 
or  to  the  propagation  model.  An  explanation  has  not  been 
determined  for  the  (J- shape  nature  of  the  wavenumber  spectrum 
plots  and  thus  slightly  different  curves  might  rs5ult  if 
this  is  found  to  bias  the  null  spacing.  Also  there  is  a 
bias  associated  with  the  method  of  manually  estimating  the 
null  spacing.  This  might  be  resolved  with  a  computer  proce- 
dure to  locate  the  nulls.  The  experimental  curves  actually 
agree  well  with   theory   under  the   assuaption   of  a   flat. 
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■^ot ally-reflecting  surface  and  a  fully-absorbing  bottom. 
The  premise  that  the  theory  might  be  approximately  accurate 
ip-  realistic  cases  appears  valid.  This  overall  generally 
smooth  fit   supports   the    validity    of  the   wr. 

This  has  been  an  initial  investigation  of  the  WT,  and 
thus  only  a  broadbrush  study  was  possible.  These  prelimi- 
nary results  indicate  that  additional  analysis  might  be  ben- 
eficial. The  response  of  the  WT  to  input  from  different 
propagation  models  would  deteriine  if  results  from  different 
models  could  be  compared.  Detailed  investigation  of  the  spe- 
cific effects  of  various  parameters  would  lead  to  a  quanti- 
tative analysis  of  the  sensitivity  of  the  WT.  An  explanation 
for  the  amomalous  U-shaped  speotrum  and  rhe  development  of 
ar.  algorithm  which  can  accurately  determine  null  locations 
would  be  necessary  before  the  »T  could  be  used  practically. 
Ar.  extremely  important  investigation  would  be  the  comparison 
of  analyzed  actual  signals  to  predictions  from  the  WT  to 
determine  applicability,  accuracy,  cop.sistency  and  any 
appropriate   biases. 

The  intent  of  this  research  was  to  determine  the  valid- 
ity of  the  HT  by  development  of  source  depth  determination 
curves        and        to        evaluate      the        oceanic        and        acoustic 
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sensitivities  of  these  curves.  The  results  obtained  indicate 
that  in  general  the  HT  appears  to  be  valid  wi*-h  respect  zo 
-he  Lloyd's  Mirror  effect.  The  parametric  variations  which 
were  investigated  show  that  ocaanic  variability  is  extremely 
important,  but  that  the  WT  is  flexible  for  other  parameters, 
such  as  transmission  geometry.  Ths  WT  appears  to  bs  viable 
for  operational  applications.  The  additional  investigations 
suggested  above  would  determine  the  validity  of  this  suppo- 
sition. It  is  conceivable  that  the  WT  could  be  developed 
into   an   operational   prediction   product. 
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